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Abstract. In a recent paper [J. Math. Phys. 47 082303 (2006)], Quantum Energy Inequali- 
ties were used to place simple geometrical bounds on the energy densities of quantum fields in 
Minkowskian spacetime regions. Here, we refine this analysis for massive fields, obtaining more 
stringent bounds which decay exponentially in the mass. At the technical level this involves the 
determination of the asymptotic behaviour of the lowest eigenvalue of a family of polyharmonic 
[>C . differential equations, a result which may be of independent interest. We compare our resulting 

' bounds with the known energy density of the ground state on a cylinder spacetime. In addition, 

. we generalise some of our technical results to general L^-spaces and draw comparisons with a 

' similar result in the literature. 

X) 
D 

^ ■ 1. Introduction 



Quantum Energy Inequalities (QEIs) quantify the extent to which a quantum field can vio- 
^ I late the energy conditions of classical general relativity. For example, the real scalar field in 
d-dimensional Minkowski space admits physically reasonable stated for which the expected 
O ■ energy density is negative, thus violating the Weak Energy Condition. Moreover, the energy 
^ ! density at any given point is unbounded from below. However, the theory also satisfies a QEI 
I bound [fT0l[T2l which asserts that averages of the (normal ordered) energy density along an iner- 
O ■ tial curve 7, parameterized by proper time and with velocity m", obey 



9": [ dr {u^u':Tat:{l{T)))^g{rf > -K, T ^u''\g{u)\'Q,{u/m) (1) 

^ ■ J Jm 

> ■ 

' for all physically reasonable states ip, where g is any smooth, compactly supported real-valued 
function, ^(m) = J e^^'^'^g^r) dr is its Fourier transform, 

><: Q.{x) = ^^j\yy\y^-ir^'y^ (2) 

and the constant Kd = Ad-2/ (2(i(27r)'^~^), with Ak the area of the unit fc-sphere. Similar QEIs 
are known for a variety of free field theories in flat and curved spacetimes, and also for positive 
energy conformal field theories in two-dimensional Minkowski space. We refer the reader to the 
reviews [|9l|22l and [fTTI for references and applications of QEIs. 

For many purposes it is convenient to have a simpler form for the QEIs. One way of doing this 
was developed recently in fillil : estimating {u^u'':Tab-i'j{T)))^ by its supremum over an open 
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^By 'physically reasonable' , we mean Hadamard states, which have smooth normal-ordered two-point functions. 
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interval / C M, we have 

r ^u^\g(u)\^QJu/m) 
sup(«V:T„,:(7(r)))^ > -i^ /^ - Jf , ^ ^ (3) 
re/ j ar5((r)^ 

for all real-valued (7 compactly supported in /. As the left-hand side is independent of g, we 
are free to optimize the inequality over the class of permissible g. For massless scalar fields in 
even-dimensional Minkowksi space, this may be converted to an eigenvalue problem, leading to 
the bound 

snp{u^u':T,MT)))^ > (4) 

re/ 

where tq is the length of / and A„ is the smallest eigenvalue of the polyharmonic equation 

= m) (5) 

on (—1, 1) subject to boundary conditions ^/'(±1) = = ■ ■ ■ = ?/>("~^)(±l) = 0. The same 

result applies to fields of mass m > (with the minor change that an overall factor of 6/5 must 
be inserted on the right-hand side of (HI) in d = 2 dimensions). However, as was noted in ifTTI . 
this is a rather weak estimate if mro ^ 1: fixing any g, it is easy to show that the right-hand side 
of ^ decays faster than any inverse polynomial as mro 00. Thus we know that bounds of the 
type 

sup(«V:T„,:(7(r)))^ > (6) 

re/ V^'^Q) 

exist for suitable constants Cd,n and any integer n > d/2. This information in itself is not 
particularly useful, unless supplemented with a discussion of how the constants Cd,n grow with 
n. The purpose of the present paper is to investigate this question. In fact we will find that 
(nearly) exponential decay can be obtained; our main result is that the above bound is satisfied 
with 

22"+ir(n + l)r(2n + 1/2) 
T{n + 1/2) 

for n > d/2. Here K'^ = Kd except for d = 2, where K2 = 6K2/5. 
We are also free to optimize over n. This procedure leads to a bound 



(~'d,n = J^d T./,, , , /ON (') 



sup(MV:r„b:(7(r)))^ > -Q(m, tq) (8) 

re/ 

where 

Q(m, To) ~ K'^2''+^^m\mTof'^e~'^^''^ (9) 

for rriTQ ^ 1. 

One of the key steps in our argument is to show that the eigenvalue A„ obeys 

2r(n+l)r(2n + l/2) ^ 2V{n + l)T{2n + 1/2) 

r(n + l/2)/2„ ^^"'^ r(n + l/2) ^'""^ 

for n > 1, where i?„ = 3^2(1/2, 1, — n; 1 — 2n, n + 1; 1) is given in terms of the hypergeometric 
function 3F2 [1]; moreover, _R„ — > 1 as n — 00, so both bounds are asymptotic to A„ in this 
limit, as is the simpler formula \/2{2n)\ (which we do not claim to be a bound). We were not 
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able to locate this fact, which may be of independent interest, in the literature. Some analogous 
results for the same operator, but with different boundary conditions, are known; see Sec. [6l 

Clearly, the results presented here represent a considerable improvement on and permit 
many of the results of [11] to be strengthened. The main thrust of [TT] is that the Minkowski 
space results just mentioned also apply in curved spacetimes, provided the segment of the inertial 
curve 7 parameterized by / may be contained in a 'sufficiently large' region in which the metric 
is Minkowskian. An example will be given in Sec. [51 

The paper has the following structure: in Sec. [2l we reduce our problem to the eigenvalue 
problem ([5]); estimates for the eigenvalues A„ are obtained in Sec.[3]and the optimisation over n 
mentioned above is performed in Sec. IH Sec. [5]contains an example in which our bound may be 
compared against a known value for the expectation value of the energy density. In this particular 
instance, the ratio of the actual energy density to our bound tends to zero exponentially as the 
mass increases so we cannot conclude that our bounds are asymptotically sharp. (Equally, we 
cannot rule out this possibility.) They nonetheless represent a distinct improvement on earlier 
results. In Sec. [6] we consider the underlying reasons for the success of the strategy employed 
in the previous sections. A mixture of theoretical and numerical evidence suggests that the 
sequence of solution operators to dS]) may be 'asymptotically rank 1'; a phenomenon which 
has been established elsewhere for solution operators to the same differential equation but with 
different boundary conditions. Finally, in Sec. |7] we show how our analysis may be extended to 
determine the L^-operator norms on the solution operator to ([5]); this is of independent interest 
and allows some comparison between our main results (in an L^-context) and a result obtained 
in |l6l (which is related to the version of our problem). 



As initial preparation, we notice that the integrand in (|2l) is bounded from above by y'^~^ 
if d > 2, and hence Qd{x) < 1 for x > 1. In the case = 2, it may be shown [|TT]| that 
Q2{x) < 6/5 on this domain. Accordingly, ^ implies 



for all real-valued g e C^{I), where K'^ = Ka for d > 2, K'^ = 6K2/5. Since the right- 
hand side is invariant under translations in g, we may assume without loss of generality that 
/ = (-ro/2, To/2). Writing g^r) = /i(2r/ro), where h e 1), we obtain 



2. Reduction to an eigenvalue problem 




(11) 



(12) 



for all real- valued h G ( 
Defining, for any such h. 



1, 1), where x = mro/2. 



/i dt\h{t)\^ ' 
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our problem is now to estimate from above the function 

X ^ inf Hd,h{x) 

h 

for X ^ 1/2. Writing {Df){t) = id// dt, observe that for any n > d/2, and x > 



X 

„d-2n / fCT, \|2 



-oo 27r 
•1 

„d-2n 



= a;"y dt|(D"/i)(t)|' 

using the monotone decrease of y'^^^n ]]^+^ ^j^g j^^^^ ^j^^^ is even because h is real- 

valued, Parseval's identity, and the fact that h is supported on (—1, 1). Introducing the usual L^- 
inner product (-, ■) on (—1,1) (by convention, this is linear in the second slot) and its associated 
norm || ■ ||, we can write the last expression in the form 

M < dt|P"/^)(t)P _ {D-h,D-h) _ ., Jh,D^-h) 

\\h\\^ ~ {h,h) ~ {h,h) 

(using the symmetry of on 1, 1)) and minimise the right-hand side over h (excluding 

the identically zero function). By Theorem X.23 in [|20ll . the infimum is the lowest element of the 
spectrum of the Friedrichs extension Aof D^" on (—1,1), whose domain is the intersection of 
Sobolev spaces [2] D{A) = W^'^i-l, 1) n Vr2"'2(-1, 1) [see, e.g.. Sec. II.B in [13]]. Moreover, 
the operator A has compact resolvent, by a straightforward modification of the proof of Theorem 
XIII.73 in En . so A has purely discrete spectrum. Using elliptic regularity, the eigenvectors 
of A are smooth solutions to Q, and since they belong to Wq''^{—1, 1) they obey the boundary 
conditions V^(±l) = V^'(±l) = ■ ■ ■ = V^("-i)(±l) = 0. 
To summarise: we have established that 

inf/frf,;,(x) <x"A„, 

h 

where A„ is the minimal eigenvalue of ([5]) subject to the boundary conditions just mentioned. 

3. Estimates for the Minimal Eigenvalue 

Let G : [—1, 1]^ M denote Green's function for ([5]), and T denote the associated solution 
operator 

{Tfm = l'^G{t,s)f{s)ds. 

Since this is the inverse to the original problem, we seek the maximal eigenvalue, or spectral 
radius, of T, which we shall denote r(T). Numerical investigation suggests that, for large n, the 
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eigenfunction associated with the maximal eigenvector is well approximated by f{t) = (1 — t^)"; 
this observation leads us to rigorous bounds via the following fact. 

Lemma. If a and b are positive constants such that af{t) < {Tf){t) < bf{t)for all t G [—1, 1], 
then a < r{T) < b. 



This is part of the general theory of order-preserving operators (see, for example, [1161 Lemmas 
9.1, 9.4]), but for the reader's convenience we include the short proof. 

Proof of Lemma. Since < af{t) < (T/)(t), we can square and integrate to give ||a/|| < 
and hence a < \\T\\, where ||T|| denotes the operator norm of T acting on 1, 1). Since T is 
self-adjoint, ||T|| = r(T), so we have a < r(T). Note that the operator norm inequality is true 
for any norm, not just forp = 2; we exploit this fact in Sec.|71 

The image under T of any function is 2n — 1 times differentiable and has derivatives of 
order up to n — 1 equal to zero at both endpoints; in particular, it can be written as p(t) = 
(1 — t^Yqif) = f{t)q(t), where q E C[—l, 1]. The Banach space X of all such functions, with 
norm 

II II I^WI 

lbll/= sup — — 

te{-i,i) JW 

is therefore T-invariant and contains all of the eigenfunctions of T. In particular, the spectral 
radius of T as an operator on X is the same as its spectral radius as an operator on 1, 1). 
For any p G X, we have by definition 

-\\p\\ff{t)<p{t)<\\p\\ff{t). 

Since G is non-negative [4], we can apply T to this inequality to give 

-MfiTfm < {Tp){t) < \\p\\f{Tf)(t). 

Since {Tf){t) < bf{t), we have 

-b\\p\\ff{t)<{Tp){t)<b\\p\\ffit) 

which is to say that ||Tp|| j < This shows that the operator norm of T on X is no larger 

than b, and hence that r(T) < b. □ 

Note that the lower bound on r (T) does not depend on the positivity of the Green function G. 

To find suitable constants, we shall find an exact formula for Tf. Although there is an explicit 
formula for the Green function its use would involve some integrals which are not obviously 
tractable; we therefore exploit the fact that / and Tf are both polynomials of known degree, to 
reduce the differential equation to a finite system of linear equations. 

Since / is a polynomial of degree 2n, all of its (2n)th order integrals are polynomials of 
degree An, exactly one of which satisfies the boundary conditions, or equivalently has a factor of 
(1 — t^)". Moreover, the differential operator and boundary conditions commute with reflection 
in the origin, so the same is true of T; since / is even, Tf is also even. In view of the factor 
(1 — t^)", it is convenient to write (T/)(t) = (1 — t^)"P(l — t^), where P is a polynomial of 
degree n, say P{z) = Ylk=o ^kz'^ ■ The bounds for the spectral radius are then 

min Pil-t^) < riT) < max Pil-t^). 

te[-iA] te[-i,i] 



SIMON P. EVESON AND CHRISTOPHER J. FEWSTER 



To determine P, we must solve the equation 



J2n " 
fc=0 



2\n+k 



[l-t 



2\n 



We appproach this simply by expanding the powers of 1 — using the binomial theorem, differ- 
entiating, and equating coefficients. 



:-ir 



dt2'^ 



n + k 



k=0 r=n 
n k 



n + k 
r 

(2r)! 



2r 



2\n 



r J [2(r — n)]\ 



2(r-n) 



2\n 



fc=0 j=0 

(substituting r = n + j) 



n + k\ [2{n + 



n + j 



j=0 



.k=j 



n + k\ [2{n + j)]\ 



n+j 



(2j)! 



2\n 



'l-t 



2\n 



We can now equate coefficients of on each side to give the equations 

'n + k\[2{n + 
n + j 



k=j 



(2j)! 



(0 < J < 



In matrix form, these read ABa = (3, where (3 is the vector ^(— 1)-' (") j , A is the diag- 
onal matrix diag ^((— 1)-' [2(r7, + j)]!/(2j)!)J^pj , and B is the matrix of binomial coefficients 
f ("+*^)) ; here ("i'') is understood to be zero if n + j > n + k. There is a simple for- 
mula for the inverse of this matrix, following from the identity Er=p(-l)^^''C) (D = ^pq- 
B-^ = [(-1)^+'= (;;+^) ) " .We therefore have a = B'^A'^p, or explicitly 



k=j 



n + j J [2{n + k)]\ \kj' 



We simplify this expression following the procedure in [[T9l §3.3]. Letting r = k — j,we have 



a,- 



r=0 



r2 + j + r + l\ (2j + 2r)! 



n 



The first term (r = 0) is 



n + j J {2n + 2j + 2r)\ \j + r 
{2j)\ [n 



[2{n + j)]\\j 
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and the ratio of term r + 1 to term r is 

(r + j + l/2)(r + j - n) 
(r + j + n + l/2)(r + l)' 

(note that this formula gives if r = n — j). We can therefore identify the sum as a hypergeo- 
metric function 



«i = ro/!, '-M, ( . ) 2^1 (j + 1/2, j-n;j+n + 1/2; 1). 



(2j)! fn 

Gauss's identity |[T9l §3.5] states that, provided Re(c — a — b) > 0, 

, , T(c)T(c-a-b) 

,F,{a,b;c;l)= ) ' \ -(. 

r(c — a)r(c — b) 

We apply to this to give 

T{j + n + l/2)T{2n-j] 



2 



Fi(j + l/2,j-n;j+n + l/2;l) 



r(n)r(2n + 1/2) 
and can therefore conclude that 

nr(l/2 + n + j)T{2n - j)V{l + 2j) 
r(l/2 + 2n)V{2n + 1 + 2j)r(l + ])V{n + 1 - j) ' 

Since each aj is positive, P(l — t^) attains its minimum and maximum over [—1,1] at ±1 and 0, 
respectively, and we have 

n 

ao = P(0)<r(T)<P(l) = 5^a,. 

i=o 

From the point of view of the application, the lower bound is the more important one (and, as 
already mentioned, does not depend on the positivity of the Green function). Before we consider 
this, though, we shall show that the ratio of the upper and lower bounds tends to 1 as n ^ oo, so 
both bounds are in fact asymptotically equal to r(T) as n ^ oo. The ratio is 



1 " 

Rn = ^ 



and we can use the same technique as before to identify this sum as a hypergeometric function. 
The ratio of two successive terms is given by 

(2j + l){n- j) (j + l/2)(j - n){j + 1) 



so 



aj 2{2n - 1 - j)(n + j + 1) (j + 1 - 2n){3 + n + l)(j + 1) 



1 " 

- 5^ «i = 3^^2(1/2, 1, -n; 1 - 2n, n + 1; 1). 



an 

" i=0 
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We can now calculate the limit; our strategy is influenced by an unpublished calculation of T.H. 
Koornwinder, also employed in ifTSl . We first expand in Pochhammer symbols to give 



R 



l/2)k{l)k{-n)k 



where the sum terminates at n because {—n)k = for k > n. The (1)^ term in the numerator 
cancels with the k\ term in the denominator, and the other symbols can be expanded to give 



k 
k=Q 



where 



^ [(-n)(l - n) ... (A; - 1 - n)][(l/2)(3/2) ...{k- 1/2)] 

[(1 - 2n)(2 - 2n) ... (A; - 2n)][{n + l)(n + 2) . . . (n + k)] 
1 [n(n A; + 3.5... (2fc-l)] 

~ ¥ [{2n - l)(2n - 2) . . . (2n - k)][{n + l)(n + 2) . . . (n + k)] 

< — =: dk 
-2k 

The last step is true because each term on the numerator is no larger than the corresponding term 
on the denominator; specifically, 

n — r + 1 < 2n — r; 2r — 1 < n + r (1<t< n). 

We also have c„ q = 1 and, for each A; > 0, c„ ^ ^ as n — > oo (because Cn^k is a rational function 
of n whose denominator has degree k greater than its numerator). We now have < Cn,k < dk, 
Yl'kLi dk < oo and, for each k, Cn,k ^fco as n ^ oo. It follows from Tannery's theorem (i.e., 
the Dominated Convergence Theorem on a measure space of countably many atoms of mass 1) 
that 

n oo 
k=0 k=0 

We now know that 

r(n + l/2) 



2r(n + l)r(2n + 1/2) 



is a lower bound for r(T), asymptotically equal to r(T) as n ^ oo. It follows from Stirling's 
formula that 

r(T) ~ 

^ ^ V2(2n)! 

as n — » oo; but this is greater than ao for large n, so we cannot conclude that this is a lower 
bound for r(T). 
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4. Optimisation over n 
We know from the calculations in the previous section that 

... / X r, o„2r(n + l)r(2n + 1/2) 

We now optimise over n for fixed x, allowing, for the moment, non-integer values of n. The 
logarithmic derivative is 

F\n) = -2 log(a;) + ^(n + 1) + 2^(2r2 + 1/2) - ^(n + 1/2) 

where \1/ is the digamma function. We seek a critical point of F, so wish to solve the equation 

log(x) = + 1) + ^(2n + 1/2) - ^^{n + 1/2). 

The right-hand side has an asymptotic expansion 

log(2n) + + 0(1/^2) 
An 

(a straightforward calculation from ffTl Equation 6.3.18]) so we have 

X = 2nexp(l/(4n)) exp(0(l/n^)) 
= 2n{l + l/(4n) + 0{l/n'^)){l + 0{l/n^)) 
= 2n + l/2 + 0{l/n). 

It is apparent that n ~ x/2 as n ^ oo, so 0{l/n) = 0{l/x) and we can solve the equation to 
give 

n = x/2-l/A + 0{l/x). (13) 
Next, we have (by definition of the polygamma function and using [|Tl Equation 6.4.12]) 

F"{n) = i!i{n + 1) + 4V^i(2n + 1/2) - ^i(n + 1/2) = 2/n + 0{l/n^) 

At the critical point, (fTSl) is true; multiplying by A/{nx) (and remembering 0{l/n) = 0{l/x)) 
yields 

2 4 

- = - + 0(l/x2) 
n X 

so, at the critical point, F"{n) = A/x + 0(l/x^). We can therefore expand F about its critical 
point no {F is analytic in the right half -plane) to give 

F(no + 5) = F{n,) + 5^(1 + 0{l/x')\ + ^i^'"(^o + C) 



where C lies between and 5. We also have 



2 



F"'{n) = - + 0(1/^3) 

n^ 

(using [[ll Equation 6.4.13]) so, if we assume \5\ < 1, which implies that \C\ < 1, we have 
F"'{no + C) = 0(1/ x^) (uniformly in 6 and Q. This gives 

2P 

F{no + 6) = F{no) + + 0{l/x^). 

X 
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Now choose 6 such that \6\ < 1/2 and no + 5 G N, so exp(F(no + S)) is the bound we seek. 
Since 26"^ /x + 0(l/x^) — as x — oo, we have exp(F(no + 5)) ~ exp(F(no)) as x ^ oo. 

Finally, since moving by a distance of up to 1/2 from the critical point no has no effect apart 
from a multiplicative factor converging to 1, moving from the exact critical point uq = x/2 — 
l/4 + 0(l/x) to the approximate critical point x/2 will have no more of an effect. We thus have 

exp(F(no + 6)) ~ exp(F(a;/2)) as x ^ oo 

and we can use Stirling's formula on exp(F) to give an asymptotic formula for the bound: 

inf Hd,h{x) < exp(F(no + 5)) ~ 2^/^Tx'^^'^^'^e~'^ as x ^ oo 



5. Example: Cylinder spacetimes 

Consider the cylinder spacetime {N, r]) formed by periodically identifying four-dimensional 
Minkowski space under a translation in the 2;-direction. The ground state energy density in this 
spacetime, for the quantized minimally coupled scalar field of mass m is [|23l[T7l 

oo 2 
n=l ^ ' 

where L is the periodicity length. For mL ^ 1, the series is dominated by its first term, and so 

^^"^ ^ "(27r)3/2(mL)5/2 ^^^^ 

Now consider a timelike curve 7 : (0, To) — > given by 7(r) = (r, 0, 0, 0), which has total 
proper duration ro. If ro < L we may enclose the curve in an open globally hyperbolic subsej] 
D = /+(7(0)) n /^(7(ro)) of the spacetime. Quantum field theory in D is indistinguishable 
from quantum field theory in its isometric image in the covering Minkowski space. (See |[5l [HI 
for a full presentation of this idea.) Accordingly, we may use Minkowski space QEIs to constrain 
the energy density in D. Applying our results, we obtain an a priori bound {Tu) > —Q{m, tq). 
The best constraint is obtained for tq = L; for mL ^ 1, this is asymtotically 

167r2 

Note that although our a priori bound is consistent with the known value of the energy density, 
it does not adhere to the same exponential law. The same would be true for the ground state on 
other toroidal quotients of Minkowski space. There are three possible explanations for this. The 
first of these is that the weaker bound might be needed to accommodate states other than the 
ground state. Second, it may be that the estimate made in Sec. [2l are too weak (after that point, 
all our estimates are asymptotically sharp). Thirdly, it may also indicate that our starting form 
for the quantum energy inequality, which is known not to be optimal, becomes progressively less 
sharp at large mass. 



Q{m,L) . (16) 



^These regions were called c.e.g.h.s. regions in ifTTIl : here, we use the nomenclature of Sec. 6.6 of |[T4|| . 
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6. Remarks on the role of (1 - 1^)". 

Let T„ denote the solution operator to ([5]), as defined at tlie beginning of Sec. [3l /„(t) = (1 — 
^^)", A„ be tlie greatest eigenvalue of T„, En be the corresponding eigenspace (necessarily one- 
dimensional and spanned by a non-negative function m„ of unit norm, because of the positivity 
of the Green function; see, for example. Theorem 1 1.1(b) in [T6]), and P„ be the corresponding 
spectral projection P„/ = (m„, /) m„. 

The substance of Sec. [3] is that (T„/„)(t)/(A„/„(t)) tends to 1 uniformly in t as n oo. 
It is striking that there is such a simple formula which, in this asymptotic sense, behaves like 
the eigenfunction m„; indeed, fn is the simplest function (precisely, the monic polynomial of 
minimal degree) satisfying the boundary conditions. 

One possible explanation of this is that the sequence (T„) might asymptotically have rank 1, 
so ||T„ — AnPn||/||T„|| — i> as 72 ^ oo. If this were the case then, for any function /, T„/ would 
approach the eigenspace En- Since /„ is a multiple of T„l, its behaviour would be much less 
surprising in this context. 

There is some numerical evidence that this is indeed the case: — A„P„||/||T„ || is the ratio 
of the second-largest eigenvalue to the largest eigenvalue, which can be calculated numerically; 
it appears to behave like c/n, where c is a constant in the region of 1/4. The numerical method 
starts by calculating Green's function for Equation ([5]) explicitly for any particular n; numerical 
quadrature schemes can then be employed to calculate the eigenvalues to any required degree of 
precision. In the table below, Ai and A2 are the first two eigenvalues of T„, the solution operator 
to ©. The fourth column is converging to 1, illustrating the asymptotic expressions for the 
dominant eigenvalue derived in Sec. |3l The final column appears to be converging to a limit in 
the region of 1/4, as mentioned above. 



n 


Ai 




A2 




v^(2r2)!Ai 


nX2/Xi 


5 


2.01975 X 10- 


-7 


1.04991 X 10- 


-8 


1.03652 


0.259911 


10 


2.96037 X 10- 


-19 


7.56909 X 10" 


-21 


1.01856 


0.255681 


15 


2.69890 X 10" 


-33 


4.56884 X 10^ 


-35 


1.01242 


0.253928 


19 


1.36524 X 10" 


-45 


1.81896 X 10^ 


-47 


1.00982 


0.253144 


20 


8.74731 X 10" 


-49 


1.10651 X 10" 


-50 


1.00933 


0.252995 



This conjectural asymptotic rank 1 behaviour is known to hold for a closely related prob- 
lem, in which the boundary conditions in (|5]) are changed to ip^^\—l) = (0 < j < n — 1), 
tl>^^\+l) = 0(n < j < 2n—l). The solution operator for this problem is related to the Riemann- 
Liouville fractional integration operator, and this yields an asymptotically correct upper bound 
of (ri!)^/2^""^ for the minimal eigenvalue of the differential operator [24J. Similar bounds were 
found independently in [15]; for tighter bounds, see [O. The asymptotic rank 1 property of 
the solution operators is an immediate consequence of results in [7] and (81 on iterated Volterra 
convolution operators. 

The same property can be seen to hold in another example. If the boundary conditions for 
Equation © are changed to V'^^H-l) = (0 < j < 2n - 2, j even) and ip'^^\+l) = 
(1 < j < 2n — 1, j odd) and T„ represent the solution operator, then it is easy to see that 
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Tn = T". The leading eigenvalue of Ti is simple, and it follows from the spectral theorem that 
T"/||r"|| is asymptotically equal to the associated spectral projection. 



7. Results For Other Lp spaces 
The solution operator to Q 

{TJ){t) = j' Gn{t,s)f{s) ds 

can be thought of as an operator on any of the spaces 1, 1) (1 < p < oo). Denote by ||T„||p p 
the operator norm of r„ acting on 1, 1). In Sec. [3l we found asymptotically correct upper 
and lower bounds for ||r„||2,2; denote these by 



v^(2n)! - " " ' - V2(2n)! 

where a„ and hn both tend to 1 as n ^ oo. Moreover, as remarked in the proof of the Lemma in 
Sec. |3l the lower bound is valid for all norms, so in fact we can write 



\/2(27z)! 



< IIT I 



for all -p e [1, oo]. It is easy to calculate exactly ||T„||oo,oo, as follows. By definition, T„l = gn, 
where gl^'"^ = 1 and gl^\±l) = (0 < j < n - 1); it follows that = (1 - t2)"/(2n)!. 
Since ||(7„||oo = l/(2r2)! and ||l||oo = 1, we have a lower bound HT^Hocoo > l/{2n)\. Moreover, 
for any / G 1) and t G [-1, 1], 

liTnfmi < ll/lloo Gn{t, S) ds = || / 1| oo (T„ 1) (t) = ^(1 " t^)"||/||oo 

(using the non-negativity of Gn)- Taking a supremum over t G [—1,1] gives ||T„/||oo < 
||/||oo/(2r2)!, so ||T„||oo,oo > l/{2n)\. In combination with the previous inequality, this shows 
that ||T„|U,oo = l/(2n)'!. 

It now follows from the symmetry of Gn that ||T„||i i = l/{2n)\ (because T„ acting on 
L°°(-l, 1) is the adjoint of T„ acting on ^^(-1, 1)). 

Information about ||T„||p p for other values of p can now be obtained from the Riesz-Thorin 
interpolation theorem (Theorem IX. 17 in [20]). A special case of this, informally stated, is that 
if T is bounded on L^" and L^'^ and p^^ = (1 — u)pq^ + up^^ (0 < u < 1), then T is bounded 
on and \\T\\p^,p^ < ||T||i--J|T||-^^p^. Withpo = 1 andpi = 2, we have p-^ = 1 - m/2 and 

IIT II 



or, in terms of some p E [1,2], 

IITII < 52/9^ L 

U\\p,P-0^ (2n)!2V</ 
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where p ^ + q ^ = 1. Similarly, with po = 2 and pi = oo, we have p^^ = (1 — u)/2 and 



IT II < ft^"" ^ ^ 



(2n)!2(i-«)/2 



or, in terms of some p G [2, oo] 



IT II < ^— 

|in||p,p^o„ (2n)!2Vp 



In general, for any p E [1, oo], we have 

IIT II < h'^l''— ^— 

||i„||p,p^o„ (2n)!2iA 

where r = max(p, q). As mentioned above, we also have, for any p, the lower bound 

^ < IIT II 

x/2(2n)! - " """'^ 

so the decay rate of ||T"||p_p is, up to a constant, independent of p: ||T„||p_p x l/[2n)\. 
Finally, we note that the identity 

inf = 2'-V! (rGM) (17) 

h£9(-l,l) 
fh=l 

is obtained in Appendix C of [6J. Although there are similarities to the p = 1 case above, there 
are also significant differences: it is an extremum over a hyperplane, as opposed to a ball, and 
there are no explicit boundary conditions. Our result permits us to prove the related bound 

inf ||/.(2")|K> ||T„||^i = (2n)! (nen) 

h€&{-l,l) 
I\h\=l 

which is weaker than the result of [6] (for r = 2n) by the geometric factor of 2^"^^ The 
origin of this is likely to be the absence of boundary conditions in that result, as it seems that 
one may approach the bound by nonnegative h (so the difference between the integral and 
norm is inessential). Boundary conditions enter because our result could equally be stated as the 
infimum as h varies over the range of T„ in L^, all elements of which obey the specific boundary 
conditions we have imposed. In fact, (fTTI) was used to obtain a result quite similar in spirit to our 
main result: namely, that 

inf sup \y'h{y)\ < iv^e^/^x^'+^Z^ -x/2 
fh=i 

for k E Nq and x > max{2A;,2}. (We have adapted the formula given in |l6l to our own 
conventions.) Apart from the difference in the extremisation domain, this could be thought of as 
an — i> L°° version of our result that 

inf ( r ^\y'^h{y)A"^ = inf aM^. h(x)/2 < Tr^/V+^/V^/^^ 
S\h?=l S\h\^=i 
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While the strategy employed in |l6l overlaps in part with ours, the key portions of the two argu- 
ments (leading to (fTTI) in L6J, or our bounds on A„) are quite different. 
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